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VAPOR INGESTION IN CENTAUR LIQUID-HYDROGEN TANK 


by Eugene P. Symons 
Lewis Research Center 

SUMMARY 

An experimental investigation of the vapor- ingestion phenomena, as it occurs in the 
Centaur liquid- hydrogen tank was conducted. The equations presented herein can be 
used to predict the anticipated fluid height above the tank outlet at which vapor ingestion 
may be expected to initiate as a function of the dimensionless grouping W/ (1 +B), 
where W is the outflow Weber number and B is the Bond number. Experiment data 
were obtained in small-scale Centaur liquid-hydrogen tanks at normal gravity and com- 
pared with predicted results. The overall agreement of analysis and experiment is 
very good. Additionally, predictions of minimum liquid hydrogen levels required to 
prevent vapor ingestion when restarting the engines in space and the quantities of liquid 
remaining in the tank at vapor ingestion during main engine firing are presented. 


INTRODUCTION 

The Lewis Research Center has been studying the phenomena that occur during 
liquid outflow from various tank systems. Earlier studies (refs. 1 and 2) examined the 
distortion of the liquid- vapor interface during the draining of a cylindrical tank in a 
weightless environment. More recent studies have examined the vapor- ingestion phe- 
nomena in cylindrical containers having both flat (ref. 3) and hemispherical (ref. 4) 
bottoms. Data from these studies were obtained both in normal gravity and in weight- 
lessness. Other investigations have examined the effectiveness of various baffles (refs. 
5 and 6) and the use of outflow rate throttling (ref. 7) as means of reducing liquid re- 
siduals (the quantity of liquid remaining in the tank when vapor enters the outlet line) in 
a low gravity environment. 

This report presents the results of an experimental study of the vapor- ingestion 
phenomena as it occurs in a scale model of the Centaur liquid- hydrogen tank. The ob- 
jective was to experimentally determine the height of the liquid free surface in the tank 



when vapor is initially ingested into the tank outlet. These data and those obtained pre- 
viously in reference 8 could then be compared with an approximate analysis (ref. 9), 
which predicts the vapor- ingestion hei^t as a function of Weber and Bond number. If 
the data show that the analysis is valid, the analysis can then be used to predict the 
depth of liquid in the Centaur tank above which a successful restart of the vehicle en- 
gines could be accomplished without the ingestion of vapor resulting in pump cavitation 
and subsequent loss in engine performance. 

In the selection of the tank model and test liquids used in this investigation, some 
attempt was made to simulate the Bond numbers and Weber numbers that would exist 
during a typical Centaur engine restart sequence. This resulted in tanks having inside 
diameters of 8, 15, and 30 centimeters and two test liquids, trichlorotrifluoroethane 
and ethanol. All tests were performed in normal gravity. 


SYMBOLS 

A 2 

A area, cm 

a radius of tank, cm 

2 

B Bond number, pga /cr 

b height of bulkhead, cm 

d/hQ dimensionless distance 

2/2 5 

F Froude number, Q /tt ga 

f dimensionless function 

/ denotes function of 

g acceleration, cm/ sec^ 

H dimensionless free- surface height, h/b 

h height of free surface above base plane, cm 

J surface curvature, cm“^ 

K constant 

A. 

K unit vector normal to surface 

p pressure, N/cm^ 

Q volumetric flow rate, cm /sec 

R radius, cm 


2 



r radial coordinate 

r unit vector in r direction 

5 length, cm 

u fluid velocity, cm/ sec 

2 2 3 

W Weber number, pQ /tt era 

Z axial coordinate 

Z unit vector in Z direction 

3 2 

^ specific surface tension, a/p, cm /sec 

6 angular coordinate 

d unit vector in angular direction 

p liquid viscosity, g/m-sec 

p liquid density, g/ cm 

ff liquid surface tension, N/cm 

cp function 

Subscripts: 

m denotes maximum 

min denotes minimum 

0 denotes point immediately above outlet 

<50 denotes point far from outlet 

1,2 denote principle local radii of curvature 


ANALYSIS 

In reference 9 Klavins presents an analysis and uses it to predict the vapor inges- 
tion height as a function of the ratio W/(l + B), where the Weber number W is a dimen- 
sionless parameter representing the ratio of inertia to capillary forces and the Bond 
number B is a dimensionless parameter representing the ratio of gravitational to 
capillary forces. This specific grouping of parameters was chosen primarily as a re- 
sult of experiment data (e. g. , refs. 3 to 5) which indicate that vapor- ingestion heists 
tend to correlate with Weber number in weightlessness and with Froude number in 
normal gravity. It can easily be observed that this grouping does reduce to Weber num- 
ber in wei^tlessness (since B = 0) and to Froude number in normal gravity (since 
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B » 1 and W/B = Froude number). This report presents a portion of Klavins’ analy- 
sis, concentrating on and somewhat amplifying that portion which would be applicable to 
the restart sequence for a Centaur vehicle (i.e, , the high Bond number regime). 

Referring to figure 1, we may write Bernoulli’s equation for a streamline on tiie 
interface between a point far from the drain (oo) and immediately above the drain (0) as 

Implicit in equation (1) is the assumption that the fluid velocity isne^igible at « 
and that the fluid is in viscid and irrotational. 

We may relate the velocity at the point above the drain Uq to the outflow rate Q 
by defining a control (or flow) area A into which the influx of fluid will be assumed to 
be normal and the velocity will be assumed to be uniform and equal to Uq, 

Thus we write 


UqA 


■/ 


dA s Q 


( 2 ) 


Further, the pressure difference p - p,. may be related to the surface curvature J. 

oo U 

Poo-P0=-‘"<'^oo- ^0> 


Substituting equations (3) and (2) into equation (1) and introducing the Weber number 

2 2 3 2 

W = pQ /tt era and the Bond number B = pga /a yield 


. 1 , 

B 2a2 ” 


» Jq> 

B 


(4) 


For convenience, the surface curvature term may be replaced by the dimensionless 
function f, defined as 


f 



a^J^- J(,) 


(5) 


Additionally, since recent studies (ref. 10) indicate that vapor- ingestion data at all ac- 
celeration levels can be correlated in terms of the group W/(l + B), we shall introduce 
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equation (5) and this group to yield 


B+1 25fVl+B/°° " 

7T a ' ' 

Now we consider equation (6) written in the following form: 


W 

B+1 


2 

2 5 


IT a 


A^(h 

oo 


ho) 


- + B 


1 + B 


( 6 ) 


(6 a) 


Making use of this form and noting the typical Bond numbers that occur during the 

0 

operation of Centaur (i. e. , 400 < Bq < 2x10 ), it can be seen that the bracketed term 
will be of the order of unity for all but extremely small values of the surface curvature 
function. 

Now, as vapor ingestion becomes imminent, hQ decreases much more rapidly than 
does h . We can use this fact by differentiating equation (6) with respect to hQ and 
then letting dh /dh„ 0. The resulting equation is 

oo U 


h 

OO 



1 M- 

A 9hQ 


-.-1 


1 /l\ 3f 

(Bf + 1) \f/aho 


(7) 


We may now combine equations (6) and (7) to eliminate h and provide a relation 

OO 

between W/ (1 + B) and the bulk liquid depth above the outlet with Bond number as a 
par ameter: 


W _ 2 A^ /l + Bf \ 

B+1 2„5 f V 1 +B/ 

As mentioned in reference 9 the results are applicable to any tank geometry subject 
to the assumptions of irrotational, inviscid, quasi-steady flow with the fluid filling the 
drain end of the tank and in the absence of vortexing. The analysis can be specified to a 
particular tank geometry throu^ the functions A and f. 

The proper choice for the control area function A depends on the liquid level in 


IIA- ( 1 \(A 9f 

A dh^ \Bf + l/\f/ahQ 
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the tank as well as the tank geometry. In reference 9 two control area functions were 
posed for the Centaur liquid-hydrogen tank, one for large liquid levels (i. e. , liquid cov- 
ering the top of the bulkhead) and another for low liquid levels (liquid level below the top 
of the bulkhead) . 

Limited experiment data (ref. 8) suggest that, for the Weber and Bond numbers that 
are to be expected during a typical Centaur engine restart sequence, vapor ingestion 
will occur at relatively low bulk liquid levels in the tank (below the top of the bulkhead). 
It is assumed that, at these liquid levels, the velocity field is essentially circumferen- 
tial and that the control area function can be approximated by two planar surfaces de- 
scribed by equation (8) and derived in the appendix. 



The expression to be utilized for the surface curvature function f is largely de- 
pendent on the contact angle of the liquid on the tank material and the Bond number. At 
large values of the Bond number, the interface is essentially flat, except very near tlie 
tank walls where the liquid contacts the wall. At low Bond numbers and for a zero con- 
tact angle, the interface is highly curved. Reference 9 examines four possible forms 
for the surface curvature function, two for B « 1 and two for B » 1, which are de- 
pendent on the initial liquid level. During a typical Centaur engine restart sequence. 
Bond numbers are so large that the interface is essentially flat, and, from the limited 
available data, bulk liquid levels in the tank at vapor ingestion are low ((h /b) < 1). In 

OO 

this case the liquid- vapor interface is in effect a flat ring with a small depression of 
depth h - h„ immediately above the outlet. The surface curvature function f corre- 

OO U Q 

sponding to these conditions is given in reference 9 as Tr‘^/8. However, an independent 
derivation outlined in the appendix yielded 

f = - (9) 

4 

Substitution of equations (8) and (9) into (6) and (7) yields the vapor- ingestion heigjit 
as a function of the parameter W/ (1 + B) . 

While it is expected that vapor ingestion in the actual Centaur liquid-hydrogen tank 
will occur at low liquid levels ((h /b) < 1) in its current operation mode, it is conceiv- 

OO 

able that future missions could be conducted at lower Bond numbers. This would, in 
turn, be reflected in greater free-surface heists at vapor ingestion and may give rise 
to instances in which (h /b) > 1. For these situations different values of the control 
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area function A and the free- surface curvature function f would be required. The re- 
lations given in reference 9 are 


A = Trh' 


eVa/ 3 y 


( 10 ) 


where 



and 


f =- 


( 11 ) 


The derivations of the preceding functions are presented in the appendix. 

Again, substitution of equations (10) and (11) into (6) and (7) would yield the vapor- 
ingestion hei^t as a function of the parameter W/ (1 + B). 


COMPUTATIONAL SEQUENCE 

The computational sequence described in reference 9 may be used to determine the 
vapor- ingestion hei^t as a function of W/(l +B). Basically, this involves choosing a 
system Bond number and a value of = hp/b. The appropriate equation for A and the 
value for f are selected based on the choice of Bond number and Hq. Using these val- 
ues, A can be determined from equation (8) or (10). (In some cases considered in 
ref. 9, f is a function of Hq-, however, for ttie range of parameters in this study it is 
constant. ) Next, the derivatives necessary in the computation are approximated. (Ac- 
tual derivatives could be calculated, but in li^t of the other approximations the method 
used in ref. 9 was employed. ) 
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^»,i^[A(Ho+0.01)-A(Ho)] 
and, since f is a constant for our studies, 



The values for the derivatives along with those for A(Hq) and f are then substi- 
tuted into equation (7) to obtain h - h^, which also prescribes a value for h /b since 
hQ/b was selected at the outset. The values of h^ - h^, AOiq), and f are then substi- 
tuted into equation (6) to obtain the parameter W/ (1 + B) . In this fashion a plot showing, 
the relation of h /b as a function of W/ (1 + B) at a point of vapor ingestion can be 

C30 

constructed. Such a plot is shown in figure 2 for B = 1000. The parametric spread for 
other Bond numbers in the range of interest (i. e. , 400 to approximately 2x10 ) is ne^i- 
gible, as was expected from an inspection of equation (6a), and thus the curve for 
B = 1000 could be used for all B > 400. 

Note in figure 2 that the predicted vapor- ingestion heig^it when using the planer 
area approach (circumferential flow around the bulkhead, eqs. (8) and (9)) differs from 
that calculated when using the spherical control area (eqs. (10) and (11)). If the region 
of interest is (h /b) <1, it is anticipated that the planer area approach is valid, while 

OO 

for (h^/b) > 1, the spherical area approach should apply. However, at values of 
(h^/b) w 1, it is not clear which form would be appropriate, and experimental verifica- 
tion would be necessary. 


APPARATUS AND PROCEDURE 
Experiment Tanks 

The experiment investigation used t ee cylindrical tanks. The bottoms and outlets 
of these tanks were scaled models of thi !entaur liquid-hydrogen tank, and the internal 
tank diameters were 8, 15, and 30 cent 3ters. The tanks were fabricated from cast 
acrylic plastic that had been polished to provide optical clarity. A flat disk having a 
diameter of 1 tank radius was positioned 1 tank radius below the pres sur ant gas inlet 
port to prevent the direct impingement of the incoming gas onto the liquid- vapor inter- 
face during draining. The tank was sealed at the top and bottom with an O-ring and was 
clamped between two stainless- steel plates by four threaded rods. (A photograph of a 
typical tank is shown in fig. 3). 
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Test Liquids 


Two liquids, trichlorotrifluoroethane and anhydrous ethanol, were employed in this 
study. The pertinent properties of these liquids are presented in table I. Both liquids 
had an essentially 0 static contact angle on the tank material in order to duplicate the 
static contact angle of the liquid hydrogen on the Centaur tank. A small amount of dye 
was added to the test liquid to improve the photographic quality. The addition of the dye 
had no measurable effect on the liquid properties as presented in table I. 


Experiment Package 

The experiment package shown in figure 4 is a self-contained unit consisting of an 
experiment tank, a pumping system, a photographic system, a digital clock, and an 
electrical system to operate the various components. Indirect illumination of the exper- 
iment tank is provided by means of a backlighting system which allowed the fluid be- 
havior to be recorded by a high-speed, 16-millimeter, motion- picture camera. The 
pumping system comprises an air reservoir, the experiment tank, a solenoid valve, 
and a liquid collection tank. The volume of the air reservoir is much greater than the 
largest volume of liquid drained from the tank during outflow, so the supply pressure 
remains essentially constant. Elapsed time during the tests is observed by reading a 
digital clock having an accuracy of ±0. 01 second. A scale positioned beside the tank 
provides an indication of the initial liquid hei^t as well as the interface hei^t as a 
function of time during draining. All the electrical components are operated through the 
control box and receive their power from rechargable nickel- cadmium cells. 


Test Procedure 

Before the flow components were assembled, the tank and all flow lines were 
cleaned in an ultrasonic cleaner to assure that the properties of the test liquids would 
not be affected by the presence of contaminants. The parts were then rinsed with dis- 
tilled water and dried in a filtered air dryer. All cleaning was done in the Lewis Zero 
Gravity Class 10 000 Clean Room (fig. 5). All parts were then assembled and mounted 
in the experiment package. Liquid was then added to the tank and all flow lines. The 
system was activated several times to purge all trapped air from the lines and checked 
for leaks. Draining was accomplished by pressurizing the air reservoir to a prede- 
termined pressure (dependent on desired flow rate) . A solenoid valve located down- 
stream of die tank outlet initiated the draining. 
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Before a test, the tank was filled to the desired level (generally about 1 tank radius 
above the elliptical bulkhead), the air reservoir was pressurized, and the camera was 
loaded. The tests were then initiated by remotely activating the electrical control box 
from a control room of the Zero Gravity Facility. (See fig. 6.) 


RESULTS AND DISCUSSION 
Comparison of Experiment Data with Analysis 

A series of draining tests were performed with thi ee scale models of the Centaur 
liquid-hydrogen tank (l/38, l/20, and l/lO scale) using two test liquids (ethanol and 
trichlorotrifluoroethane) over a range of outflow rates in a normal gravity environment. 
System Bond numbers ranged from 13.31 to 7805 and outflow Weber numbers from 
0.0062 to 35.12. The results obtained are compared with the analytical predictions in 
figure 7. Below nondimen sional vapor- ingestion heights of about 1.1, the data agree 
very well with the curve generated when utilizing the planer control area; for larger 
nondimen sional vapor- ingestion heights, the spherical control area approach more ac- 
curately predicts the draining behavior. The relatively good agreement between analy- 
sis and the experiment tends to support the assumptions previously made in determining 
the control area function. A summary of the data is presented in table II. 

In figure 8 data obtained in reference 8 are compared with the analytical predictions. 
In that study two Centaur scale models of the liquidr- hydrogen tank (1/3.67 and 1/.38 
scale) were used in conjunction with two test liquids (ethanol and water) for draining in 
a normal gravity environment. Again, predicted vapor- ingestion heights are consis- 
tently in agreement with those observed experimentally. As was the case for the exper- 
iment data taken in the present study, the data of reference 8 compare very well with 
the planar area approach for non dimen sional heights less than about 1.1. Also plotted 
in figure 8 is a single data point for a full-scale Centaur liquid- hydro gen tank which was 
obtained from a normal- gravity engine firing. 


Implication of Results 

By examining a typical Centaur engine restart sequence and calculating the system 
Bond and Weber numbers, the analytical curves could be used to determine if the vehi- 
cle's engines could be successfully restarted without ingesting vapor into the outlet line 
Ingestion of vapor could result in pump cavitation and subsequent loss of engine per- 
formance. 
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A typical engine restart sequence consists of an initial settling of the liquid propel- 
lant in a high-Bond-number environment followed by propellant outflow from the tank 
during feed- line chill down, engine cool down, and finally main- engine start, A sum- 
mary of the Bond numbers and Weber numbers associated with each phase of propellant 
outflow is given in table III, It should be noted that the calculated Bond numbers are 
based on the assumption of a vehicle mass of 18 160 kilograms (40 000 lb). 

Such an assumption would generally lead to conservative results, since during a 
restart it is quite likely that the vehicle mass would be appreciably less than this fig- 
ure because of the propellants expended during earlier engine burns. The lower initial 
vehicle mass would be reflected in hi^er system Bond numbers which, in turn, would 
give rise to lower values of vapor- ingestion height. 

The' predicted vapor- ingestion heists for the restart events presented in table III 
are plotted in figure 9, Of the three events the most critical is the engine cool-down 
phase. As shown, this results in the largest value of vapor- ingestion height of 
h /b = 0, 9 or, conversely, that an engine restart could not be accomplished without 

oo 

ingesting vapor into the outlet line if the initial liquid level were less than this value. 

It should be noted that the same curves could be used to predict the vapor- ingestion 
height (and, thus, the liquid residuals), which are to be expected during a main engine 
burn. Since the liquid- vapor interface may be expected to be flat because of the high 
Bond number, the vapor- ingestion height can be used to compute the volume of liquid 
remaining in the tank when vapor ingestion occurs. For example, if the Bond number 
was based on a vehicle mass of 18 160 kilograms (40 000 lb), the corresponding nondi- 
mensional vapor- ingestion height can be determined from figure 9 as approximately 

0.45. This, in turn, corresponds to liquid residuals of 2.45x10^ cubic centimeters 

3 4 

(8.65 ft ) or in terms of liquid mass 1. 57x10 grams ;34. 6 Ibm). 


SUMMARY OF RESULTS 

An experimental investigation of the vapor ingestion phenomena was conducted in 
l/38, 1"20, and l/lO scale models of the Centaur liquid-hydrogen tank. All tests were 
performed in a normal- gravity environment with two test liquids, trichlorotrifluoro- 
ethane and ethanol, over a range of outflow rates. Bond numbers B (based on the tank 
radius) ranged from 1331 to 7804, and Weber numbers W ranged from 0. 0062 to 35. 12. 
The data obtained from the tests were compared with an approximate analysis of the 
phenomena which yielded the vapor- ingestion hei^t as a function of the dimensionless 
grouping W/(B + 1). For dimensionless vapor- ingestion heists (h^/b) less than 1.1, 
the analysis, based on a plane-control area which assumes circumferential flow around 
the elliptical bulkhead, agreed very well with the experiment data. For larger dimen- 
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sionless vapor- ingestion heists, the analysis, based on a spherical control area com- 
pared well with the data. 

Additionally, previously obtained experiment data were compared with the analysis. 
These data, obtained using both water and ethanol in 1/3.67 and l/38 scale models of 
the Centaur liquid- hydrogen tank, again compared very favorably with the analytical 
model as did one data point obtained from an actual test firing in a full- scale liquid- 
hydrogen tank. 

Utilizing the analysis, it was determined that the Centaur's engines could be suc- 

essfully restarted after settling without the danger of ingesting vapor into the outlet, 

provided the initial liquid level in the tank was greater than 90 percent of the height of 

the bulkhead. The most critical phase of the Cantaur restart sequence was found to be 

the engine cool-down phase. Additionally, for main engine burn for the actual Centaur 

vehicle, it was calculated that vapor ingestion would not occur until the liquid volume in 

5 3 

the tank had dropped below 2.45x10 cubic centimeters (8.65 ft ). 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 29, 1977, 

491-02. 
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APPENDIX - DETERMINATION OF CONTROL AREA AND SURFACE 


CURVATURE FUNCTIONS CONTROL AREA FUNCTION 

As stated in the text, a proper choice for the control area function A depends on 
both the liquid level in the tank as well as the tank geometry. In reference 9 two 
choices were posed for the,control area function appropriate for the Centaur liquid- 
hydrogen tank; one for large liquid levels (i. e. , liquid level covering the top of the 
ellipsoidal bulkhead) and another for low liquid levels (i. e. , liquid levels below the top 
of the ellipsoidal bulkhead). 

If the initial liquid level is significantly above the top of the bulkhead, the velocity 
potential is generally approximated by a point sink, and an appropriate control area 
would be a spherical surface surrounding the tank drain. Such a control area is indi- 
cated in figure 10. Since it is necessary to know the surface area of this control sur- 
face for computational purpose, we will first develop a general expression for deter- 
mining surface areas and then specify that general expression to the Centaur liquid- 
hydrogen tank. 

We first suppose that there exists some arbitrary three-dimensional surface, say 
<p = ^(r, 0, Z) and attempt to find the projection of that surface onto the r-0 plane. 



We may then write, in general, 


dA(Z • K) = r dr d0 


(Al) 
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where A is the total surface area of the arbitrary three-dimensional surface, K is 
the unit vector normal to the three-dimensional surface, and Z is the unit vector in 
the Z-direction. Now 


K =-r^ (A2) 

I I 

where (p is ttie equation describing the three-dimensional surface. And 


0r r 80 8Z 


(A3) 


with r the unit vector in the r direction and 6 the unit vector in the 6 direction. 
Since 


Z = /(r, 0) 


then 


cp= Z- f(r, 0) 


And we can then express equation (A3) as 


8r r 80 


with 



(A4) 


(A5) 


(A6) 


(A7) 
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Then 



(A8) 


And substitution of equation (A8) into (Al) 



r dr d9 


(A9) 


Equation (A9) is true for any arbitrary surface and can now be specified to the 
geometry of the Centaur liquid-hydrogen tank in which we use a spherical control sur- 
face surrounding the outlet line and wish to determine the surface area of that surface. 

Now consider the origin of the spherical control surface to be located in the base 
plane (r, 6) of the Centaur liquid- hydro gen tank along a line drawn from the coordinate 
origin of the ellipsoidal bulkhead through the centerline of the vertical section of the 
side outlet at that point where the cylindrical tank sidewall intersects the line. (See 
fig. 10.) The radius of the spherical control surface is taken to be hg, and the equa- 
tion to describe the total spherical surface in cylindrical coordinates is 


2 

r 



(AlO) 


from which 


(All) 


After taking the appropriate derivatives of (All) and substituting into (A9), we 
finally obtain 
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(A12) 


All that remains is to establish the limits of integration in the r-d plane. 
Consider equation (A12) written in the following form: 



r dr d6 


1 - ( — 


(A13) 


Performing the integration over r yields 

►tt/2 


A = 2h; 



1 - 



d0 


(A14) 


min 


While it is conceivable that the function could be determined analytically, 

Klavins (ref. 9) approximates the radical in equation (A14) by a quadratic in 9, that is, 



(A15) 


With the constants to be evaluated from the conditions 
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2a 


u 




(A16) 



where the prime denotes differentiation with respect to 0 and 


d 


h 


0 



- 1 


If we now let 0, the constants 

min 

substituted into equation (A14), which, 


Kj^, Kg, and can be evaluated and 

after integration, gives 


r 


a = a2<;i. 




6a 



(A17) 


which is the area of the control surface for (h /b) > 1. 

CO 

If the initial liquid level is below the top of the ellipsoidal bulkhead, we would anti- 
cipate that the flow would primarily be circumferential around the bulkhead, and we 
could then use two planar surfaces for control areas. (See fig. 11.) In this case the 
problem is greatly simplified and we simply perform the integration directly at two 
equal values of 0 on either side of the outlet centerline. The appropriate integral is 
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A = 2 




dr dz 


which can be immediately integrated to yield 


(A18) 


A = 2ah^ - 



(A19) 


which is the desired area of the control surface for (h /b) < 1. 

oo 


Surface Curvature Function 

As stated in reference 9, the proper choice for a surface curvature function is 
strongly dependent on the magnitude of the acceleration environment and the liquid- solid 
contact angle. At high Bond numbers it is expected that the liquid- vapor interface will 
be essentially flat except at the tank wall where the liquid- solid contact angle will be 
maintained. At low Bond numbers and for near zero static contact angles, the interface 
is highly curved. As stated previously in the text, Klavins (ref. 9) examines four pos- 
sible forms for the surface curvature function f, two for Bq « 1 and two for 
Bq » 1. We will restrict our discussion to those at high Bond numbers. 

In general 


J = 




(A20) 


where and Rg are the principal local radii of curvature at any point. 

We consider first the case in which the initial liquid level is high in the tank (see 
fig. 12). At points far removed from the dip region, the interface is flat, that is, 

OO OO 


J « 

OO 



(A21) 
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While immediately at the dip region, we will assume that the dip is symmetric, that is, 

R- = R„ . This will not strictly be true because of the curvature of the tank walls. 

^0 2 
Klavins (ref. 9) makes the assumption that (neglecting terms of order (h - h„) and 

oo U 

higher) 


2 

R =R « E 

^0 ^0 2(h - h.) 

OO U 


(A22) 


Then 



4(^00- V 


(A23) 


And making use of equation (A23) in the expression for tbe surface curvature func- 
tion f (see eq. (5)) we obtain 


f«- (A24) 

4 

which is the expression for the approximate surface curvature function (h^/b) > 1. 

In similar fashion we may consider the case in which the liquid level in the tank is 
below the top of the ellipsoidal bulkhead ((h /b) < 1). (See fig. 1.3.) Again, at points 

oo 

far removed from the drain, the interface is relatively flat, thus 

R. = R_ w oo (A25) 


which leads to 


J 

OO 



(A26) 


In this case, immediately above the outlet at the dip region, the liquid- vapor inter- 
face is a flat ring which is depressed by an amount (h - h«) over a portion of the tank 
above the outlet. If we further assume that the variation in free- surface height is linear 
along radial lines eminating from the center of the ellipsoidal bulkhead at any fixed 
value of d we may say 
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(A27) 


GO 

^0 

Also, if we assume that the maximum dip at the free- surface occurs just above the 
center of the outlet line and further that the depression subtends an angle of approxi- 
mately 90*^, it can be shown that 


Ro 


^0 4(h^-ho> 


(A28) 


Making use of equations (A25) to (A27), we then arrive at 


J 

oo 



2 11 



^<^00 - V 
2 


(A29) 


which, when substituted into the expression for the surface curvature function (eq. (5)) 
gives 

f = - (A30) 

4 

which is the approximate surface curvature function for (h /b) <1. It should be noted 

OO 

that this approach differs from that taken by Klavens. His result for (h /b) < 1 was 

O OO 

7tV8. 
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TABLE I. - PROPERTIES OF EXPERIMENT LIQUID 


[contact angle with cast acrylic plastic in air, 0^~\ 


Liquid 

Surface 
tension 
at 20^ C, 
cr, 

N/cm 

Density 
at 20° C, 

3 

g/ eta 

Viscosity 
at 20° C, 

g/m. sec 

Specific 

surface 

tension, 

/3, 

cm^/ sec^ 

Ethanol 

22.3x10"® 

0.789 

1.2x10"^ 

28.3 

Trichloro- 

fluoroethane 

18.6x10“® 

1.58 

. 7x10“^ 

11.8 


TABLE n. - SUMMARY OF DATA 


Tank 

radius, 

a, 

cm 

Test 

liquid 

Flow 

rate, 

s' 

cm /sec 

Weber number, 
w = a® 

Bond number, 
B = ga®//3 

W/(l+B) 

Nondimensional 
vapo r- inge s tion 
height, 
h /b 

CO 

7.5 

TCTFE® 

152 

0.471 

4679 

l.OxlO""^ 

0.755 

7.5 



138.26 

.390 

4679 

8.3x10"® 

.736 

7.5 



95.59 

.186 

4679 

4.0x10"® 

.672 

4 



486 

31.8 

1331 

2.4x10“® 

1.519 




82.07 

.905 



6.8x10“^ 

.932 




261 

9.16 



6.8x10"® 

1.174 




511 

35.09 



2.6x10“® 

1.537 




9.4 

.013 



9. 7x10"® 

.493 




6.8 

.0062 

1 


4.7x10“® 

.447 

15 

Ethanol 

798.7 

.678 

7805 

8.7x10“® 

.823 




1965 

4.1 



5. .3x10"'^ 

1.1 




1517 

2.45 



3. 0x10"^ 

.856 




996 

1.05 



1.3x10""^ 

.847 




916 

.89 



1.1x10"^ 

.81 


^Trichlo rotrifluo ro e th ane . 

TABLE m. - SUMMARY OF RESTART PARAMETERS 


Event 

Bond number, 
B = ga®/p 

Weber number, 
W = Q®/ir®/3a® 

-C!b) 

Line chill down 

514 

0.0657 

-3.89 

Engine cool down 

3213 

.98 

-3.52 

Main- engine start 

1.98x10® 

6.18 

-5. 50 






Figure 1. - Centaur liquid-hydrogen tank geometry. 
Definition of terms. 
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Figure 2. - Calculated vapor ingestion height as function of W/(B + 1) for Bond 
nu mber of 1000. 
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Figure 3. - Experiment tanks. 



Figure 4. - Experiment package. 
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Nondimensional vapor ingestion height, Nondimensiona! vapor ingestion height, 



log [W/(B +1)] 

Figure 8. - Comparison of reference 8 data with analysis. 



log [w/(B +1)] 


Rgure 9. - Predicted vapor ingestion heights during Centaur engine restart. 
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Free su rface 
of liquid-\^ 


Ellipsoidal 
bulkhead — - 



Spherical control surface 

Intersection of spherical surface 
with ellipsoidal bulkhead 


t T 


^0 


Front view 


Origin of spherical control 
surface of radius hg 


Figure 10. - Approximate control area function ((fto/b) > 1). 




Figure 11. - Approximate control area function ((hoo/b) < D- 





hQ ^|o 



Figure 12 . - Approximate surface curvature function ((hoo/b) > 1 ). 
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Figure 13 . - Approximate surface curvature function ((hoo/b) < 1 ). 
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